Nonlinear propagation of light in Dirac matter 
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The nonlinear interaction between intense laser light and a quantum plasma is modeled by a 
collective Dirac equation coupled with the Maxwell equations. The model is used to study the non- 
linear propagation of relativistically intense laser light in a quantum plasma including the electron 
spin-1/2 effect. The relativistic effects due to the high-intensity laser light lead, in general, to a 
downshift of the laser frequency, similar to a classical plasma where the relativistic mass increase 
leads to self-induced transparency of laser light and other associated effects. The electron spin-1/2 
effects lead to a frequency up- or downshift of the electromagnetic (EM) wave, depending on the 
spin state of the plasma and the polarization of the EM wave. For laboratory solid density plasmas, 
the spin-1/2 effects on the propagation of light are small, but they may be significant in super-dense 
plasma in the core of white dwarf stars. We also discuss extensions of the model to include kinetic 
effects of a distribution of the electrons on the nonlinear propagation of EM waves in a quantum 
plasma. 

PACS numbers: 52.35.Mw,52.38.Hb,52.40.Db 



C/3 

'Sh! 
o 

Oh. 



> 

o\ 
rn 

O 



X 



INTRODUCTION 



The introduction of intense lasers has lead to a great 
variety of applications, including plasma based particle 
acceleration to relativistic energies [l|, Qj , and with X-ray 
free-electron lasers Q there are new possibilities to ex- 
plore dense matter on atomic and single molecule levels. 
On these length scales, of the order of a few Angstrom, 
quantum effects play an important role in the dynam- 
ics of the electrons. Using novel laser scattering tech- 
niques, quantum dispersive effects have been observed 
experimentally both in the degenerate electron gas in 
metals and in warm dense matters [J]- Hence, it is ex- 
pected that quantum mechanical effects must be taken 
into account in intense laser-solid density plasma inter- 
action experiments [5|-l7|, and in quantum free-electron 
laser systems |8l-[l(]|. Even though 7-ray lasers have not 
yet been manufactured, there have been suggestions that 
such lasers could be realized by means of annihilation 
of Bose-Einstein condensated positronium [111 . [l3 | , or by 
the excitation and nuclear spin relaxation in a lattice of 
thorium atoms [13]. This would lead to a new regime of 
intense laser-plasma interactions, where the relativistic 
quantum dynamics plays a decisive role. Intense x-ray 
and 7-ray sources exist naturally in astrophysical objects 
in the form of x-ray and 7-ray repeaters, etc. p^| - [l6| . 
In the past, the linear plasma response for relativistic 
(i.e. relativistically distributed) quantum plasmas were 
studied by deriving the longitudinal and transverse re- 
sponse functions for mildly and strongly degenerate elec- 
tron distributions [13, [3- It was noted [13] that for 
super-dense plasmas where hupe > 2meC^, there is a 
possibility of collisionless pair creation, where h is the 
Planck constant divided by 27r, ujpe the electron plasma 
frequency, rUe the electron mass, and c the speed of light 



in vacuum. The results were extended by using a Wigner 
functions approach [l9J2a, and by considering the lon- 
gitudinal response l 23l |2J] , and more general results for 
different distribution functions have also been obtained 
[25l - [27t . Relativistic quantum fluid models have recently 
been derived [2^, partially based on earlier works [20| 
of fluid-like formulations of the Dirac equation. When 
the intensity of the electromagnetic (EM) wave reaches 
a critical level (e.g. around 10"'^^ W/cm for one micron 
wavelength lasers), the relativistic electron mass increase 
and the associated nonlinearity plays a significant role 
for the propagation and dynamics of the EM wave [30] ■ 
In addition, the relativistic ponderomotive force [3l| pro- 
duces density modifications in the plasma, and the com- 
bined effects of the relativistic electron mass increase and 
relativistic ponderomotive force can lead to a modula- 
tional instability and collapse localization of EM waves 
[32, _33] . Clearly, for intense EM waves interacting with 
the plasma in the X-ray and 7-ray regimes, both rela- 
tivistic and quantum effects must be taken into account 
on an equal footing. 



In this paper, we present a nonlinear model, based 
on the Dirac equation coupled with the Maxwell equa- 
tions that are capable of treating both the relativis- 
tic (propagation and mass increase), quantum (tunnel- 
ing/diffraction) effects, and electron spin effects. The 
mathematical aspects of this system has been studied in 
the past [3J]. We will here use the basic model to in- 
vestigate the nonlinear propagation of large amplitude 
EM waves in a quantum plasma with different spin po- 
larizations. Our work has potential applications in laser- 
matter experiments [j, |35[ , quantum free-electron laser 
systems [8l4lOl]. as well as in astrophysical environments 



II. THE MATHEMATICAL MODEL 

The quantum mechanical description of the relativistic 
dynamics of an electron in an EM field is given by the 
Dirac equation 



W^p -ca-Vtp- nieC^Pip = 0, 



(1) 



where we have defined the energy and momentum oper- 
ators as 



d 
at 



and 



V 



-ihW + eA, 



(2) 



(3) 



respectively. Here, and A are the scalar and vector 
potentials, and e is the magnitude of the electron charge. 
The vector a = ax^ + ctyY + azZ, where x, y and z are 
unit vectors in the x, y and z directions, have components 
consisting of the Dirac matrices 



ttfc 



at 
<Jk 



X, y, z. 



(4) 



where the Pauli spin matrices are 



1 

1 o; ' 



and the matrix /? reads 



-i 
i 



/3 = 



I 
-I 



J-l)' (5) 



(6) 



with I being unit 2x2 matrices. 

We now wish to use the charge and current densities 
as sources for the self-consistent EM scalar and vector 
potentials for a quantum plasma. We therefore let the 
4-spinor f/; — [V'l '02 V'3 "^iC" represent an ensemble of 
electrons (T denotes the transpose of the matrix). The 
electric charge and current densities are obtained as 



-et/rtp 



Ul, — 



^e^llV'j 



and 



etp^ catp , 



(7) 



(8) 



respectively (where ip'' = [ipl 4'2 V's V'l])- The current 
density incorporates both the particle current and spin 
current. The charge and current densities obey the con- 
tinuity equation 



dt 



Vje 



0. 



(9) 



The self-consistent vector and scalar potentials are ob- 
tained from the EM wave equations 



9*2 



- c^V X (V X A) 



dt 



and 



^ dA 1 , 



MoC^je 



Pi). 



(10) 



(11) 



where /io is the magnetic vacuum permeability, Eq is the 
electric permittivity in vacuum, c = l/y^eo/Io, and pi is 
the neutralizing positive charge density of the ions. For 
immobile, singly charged ions, we have pi = eng, where 
no is the equilibrium ion number density. In our model, 
we have neglected the fact that degenerate, cold elec- 
trons are distributed uniformly in momentum space up 
to the Fermi sphere. The Fermi pressure plays an im- 
portant role in the dynamics of longitudinal electrostatic 
waves, where it contributes to the dispersion of the waves. 
For transverse electromagnetic waves, which will be our 
main interest here, the distribution of electrons play a 
minor role. The Fermi pressure is unimportant since the 
transverse electromagnetic waves are not associated with 
density perturbations. The effects on the current of par- 
ticles streaming in one direction is canceled by particles 
streaming in the opposite direction so that the net effect 
on the electromagnetic wave is negligible. For extremely 
dense plasmas, where fiiOpe is comparable to rrieC^ , the 
speeds of the electrons on the Fermi sphere become rel- 
ativistic and one can expect a frequency downshift due 
to the relativistic mass increase of these electrons. This 
effect, however, is outside the scope of our model. 



HI. CIRCULARLY POLARIZED LIGHT IN 
DIRAC MATTER 

We here consider the nonlinear propagation of light in 
Dirac matter with different spin polarizations. Solutions 
have been obtained in the past [36, 37] for single par- 
ticles in an EM field. Here we formulate the problem 
in a plasma environment, where we require the quasi- 
neutrality and current-neutrality along the propagation 
direction of the EM field. The plasma environment in- 
troduces a dimensionless quantum parameter 



flUJpf. 
^ = 2' 



(12) 



which compares the plasmonic energy Tiujpe to the elec- 
tron rest mass energy meC. Typical values are H = 10" 
for the electron number density Ue '^ 10 



30, 



in solid 

density laser-plasma experiments and H — 0.007 may 
be representative of the modern laser-high density mat- 
ter experiments |38l440| with Ue ~ 10^** m~'^. This corre- 
sponds to Wpe = 8 X 10^^ s~^ and Ae = 4 x 10^^ m for iJ = 



10- 



and LUr, 



5.4 X 10^«s-^ and Ae = 5.5 x 10"^^ m 



for H = 0.007, where Ag = c/cupe is the electron skin 



depth. On the other hand, in extremely dense plasmas 
in the core of white dwarf stares, the quantum parameter 
H may be of the order unity. For _ff > 2, it has been no- 
ticed \VT\ that there is a possibility of pair creation, and 
one then has to take into account positrons on the plasma 
dynamics. We do not consider this case here, since it is 
beyond the scope of our model. 



A. Solution of the Dirac equation 

The first step to a self-consistent picture is to solve the 
Dirac equation for a circularly polarized EM wave with 
constant amplitude. We assume a right-hand circularly 
polarized EM wave of the form 

A == ylo(xcos6' — ysinfl), where 6 = k^z — ujot, (13) 

where the frequency and wavenumber uq and feg are con- 
stants, and we assume that (f) = 0. In this case, the Dirac 
equation can be formulated into an eigenvalue problem 
(See Appendix |A| 



^0. 



fco> 



[nin + ^) - m,c^] vfi - ch{K + -f)^. 






ko. 



[h{n - y) - mec"] *2 + chiK - ^)*4 



ceAo^3 = 



- ceAo'^2 



ko. 



2 ,+mec2]*3-cM^+^)*i 



(14) 



(15) 



(16) 



and 



[n{n 



Wo 

2 
ceAo^fi 



= 0, 

) + meC^] *4 

= 0, 



ch{K 



ko. 



*-) 



(17) 



for the constant spinor components ^1-^4 where fl takes 
the role of an eigenvalue for given values of K. The 
coefficient matrix is real and symmetric, hence il is real 
and $j can be taken to be real. 

Eliminating ^1-^4 in p4)) - pT|) . we obtain the char- 
acteristic equation 



D+D_ + [uj'a 






0, 



(18) 



where 



n± 



Wo 



K± 



? 4 ? 



and we have denoted 70 = (1 + e^^§/mgC^)^/^. We note 
that D+ and D- are Klein-Gordon operators that are 
coupled in Eq. (fT8|). but which become uncoupled in the 
vacuum case wo = cko . 



B. The nonlinear plasma susceptibility and the 
dispersion relation for w^aves 

We here discuss the collective plasma response in the 
presence of the EM wave. We require quasi-neutrality 
Pe + Pi = in the plasma, which leads to 



4 

E 



IV', 






no 



* 



The current density is (see Appendix IB]) 

je =- 2ec[(^i$4 + *3*2)(xcos6'-ysin^ 

+ (^l*3-*2^4)z]- 



(20) 



(21) 



For our study, it is natural to require that the system is 
at rest in the z direction, so that jz = 0, i.e. 



vp^vl/g - ^2^i = 0. 



(22) 



The resulting current density je — — 2ec(^i^4 + 
^3^2)(xcos0 — ysin0) is right-hand circularly polarized, 
similar to the vector potential A in P^ . 

Inserting the expressions for the circularly polarized 
current je and vector potential A into the EM wave equa- 
tion (ITU)) with (j) = 0, gives 

(23) 



{lu^ - c^fco)Ao = 2^oec3(*i*4 + *3*2), 



which shows the dependence between cjq and fco- We can 
also write Eq. 



as 
,2 



2i,2 _ 
C /Cn — 



-XeWo, 



where 



Xe 



260(^1*4 + ^3*2 

eo^o^o 



(24) 



(25) 



is the electric susceptibility of the quantum plasma. 
Equation (|24l) . together with the quasi-neutrality and 
current- neutrality conditions (pil)) and ([2^ . and the 
Dirac system ([H)) - ([T71) forms a self-consistent system for 
the unknowns ^1-^4, ^, K, and ujq for given values of 
H, Aq and kg. 



C. Polar representation of the Dirac equation 

The quasi-neutrality condition (pHl) suggests that '^i- 
^4 could be represented with a polar representation. 
In addition, we wish the current-neutrality condition 
(P^ to be fulfilled. Both these conditions are fulfilled 
if we make the special choice of polar representation 



(19) *J = iV^/^)Yj with 



Yi = cos(/92 + sinipi 
Y2 ~ cos (fii + sin (p2 
Ys — cos ip2 — sin (pi 
Y4 = cosLpi — sm.ip2. 



(26) 

(27) 
(28) 
(29) 



The electron susceptibility (|25|) then takes the simple 
form 



Xe 



ecno 

2 A COs{ipi +(/32), 



and Eqs. ([T4|) ~ ([T7| become, respectively, 



W0^ 



fco^ 



- ceAoYi = 

[?i(17 - ^) - m,c2]r2 + c;i(if - ^)Y, 

- ceAnYs = 



[n{n + y) + mec^] Fg - cn{K + -j-)Yi 



2 ' 

ceAoY2 = 



and 



Wo 



[m - ^ 



)+mec2]y4 + c;i(J^-^)y2 



ceAoYi = 0, 



(30) 

(31) 
(32) 
(33) 



(34) 



which, coupled with Eq. (I24p . gives the unknowns (pi, 
ip2, r2, ii", and ujq for given values of ^Oi ^o Sind H. 
The general solution is difficult to find in terms of sim- 
ple expressions, but can be evaluated numerically with 
standard methods, e.g. Newton iterations. Some special 
choices of ipi and ip2, and the corresponding ^1-^4 are 
shown in the table below: 



(0) 



equation, we thus set Aq to zero and ^j — ^ 
P^ - P7|) . Two possible solutions of the resulting sys- 
tem are found by choosing 'i'\ 



7^) and *^°' 



= vl/(0) ^ ^ 



f^ = *o (where ^0 = 
f ^ 0, or *f = vl,„ 

and *^°^ = *f^ = ^f = 0, with K = Tfco/2 and 
n = =pajo/2 + nieC^ /fi where the upper sign corresponds 
to ^1 nonzero and the lower sign to ^2 nonzero. (Other 
solutions also exist with ^1*3 or ^4 non-zero, which 
correspond to pair states and which we, however, do 
not consider here.) Considering the first-order quan- 
tities in P^ ~ (IT71) . where we neglect first-order quanti- 
ties multiplied by each other, we find for the case with 



\l/i nonzero that ^ 



(1) 



0, 



and {—hijjQ + 2meC^)^ 
from which we have ^ 



(1) 



(1) _ 



chko^^2^ - 



(1) 



chko'9 

i(i),T,(o) 



0, 



ceAi;'-^';'' = 0, 



ccAq '^o'^o/(2wo»TieC^ 



Ticjg -I- ?ic /cq), while for the case with ^2 nonzero, we 

c;i/fco*3^' = 0, and 



instead have ^^4^' = 0, hiUQ^[^'' 



i(i) 



(1) 



{huja + 2mec2)*^^ - chko'^\' - ceA^^^o = 0, from 
which we find ^Pg ' = ceA^ ' '^ qujq / {2ujQmeC^ + huj^ — 
hc^kQ). Inserting the resulting susceptibilities Xe = 

Xe±, where Xe+ = -2ec$^°^$f V'^o'^o^o^^ and Xe- = 



-2ec*^"'*W/eou;2^y 



^,(0),^ 



,2/1(1) 



Xe± 



pe 



uJq[1 =F ^{uJq — c^kl)/2ujQmeC^ 



(35) 



into the dispersion relation (p4)) . we obtain 



Vl 


^2 


*i 


*2 


*3 


^4 


7r/2 





v^ 














7r/2 





%M) 








-7r/2 











^/na 








-^2 











v^ 



The two first lines correspond to positive energy states 
of the two spin polarizations (spin 'up' and spin 'down'), 
while the two last lines correspond to negative energy 
states, or pair states. 



D. Special cases of plasma susceptibilities 

We now consider some special cases where we can find 
simple expressions for the electron susceptibility of the 
plasma, as well as some numerical solutions of the fully 
nonlinear case. 



1. Linear and nonlinear propagation of light 

We first consider linear propagation of waves where 
eAQ/rrieC <C 1. We linearize the system by setting 



,(0) 



,(1) 



,(0)| 



,(i)| 



1(1) 



m^ = mY>+m^^\ where |^}"^| > I*]- ^|, while ^0 = ^0 
is a first order quantity. For the zeroth order Dirac 



.2^2 



z Zl 

CJg — C fc| 



W„ 



° 1 =F ?i(wo — c2/cQ)/2aJomec2 ' 



(36) 



which, after reordering of terms, takes the more trans- 
parent form 



LO, 



' c^kl-uL 



± 



nH-c^klf 



2mf,c^ujf) 



(37) 



We see that in the classical limit fi 
dispersion relation uj, 



' c kn ■ 



0, Eq. dlZl) gives the 
: for the EM waves 
in a cold electron plasma, while in the zero density limit 
-^ 0, Eq. ([57)1 yields either the vacuum EM wave 

0, or the free particle 



^pe 



dispersion relation ujq — c k\ 
equation of motion ojq ± 2mec'^ujQ/h — c^k^ = 0. The 
latter (including the pair branches) was also found from 
the electrostatic wave dispersion relation using the Klein- 
Gordon-Maxwell model |41|. 

In Fig. [T^)-c), we have displayed the solutions of the 
linear dispersion relation (j37]), and plotted the dispersion 
curves for EM waves for different values of H. Figure[T^), 
for H = 0.5, exhibits a high-frequency pair branch, the 
two EM branches shifted approximately zLH/A compared 
to the plasma frequency at /cq = 0, and a low- frequency 
branch. For H = Hcrit = 4/3-\/3 « 0.77, the pair branch 
merges with the upshifted EM branch for small wavenum- 
bers, as seen in Fig. [TJd), and for H > Hcrit the system 




ck„/(o 

pe 



ck„/(o 
pe 



FIG. 1: Dispersion curves for the linear and nonlinear prop- 
agation of light in Dirac matter for different values of H — 
fiiUpe/meC? . For the linear cases in a)-c), the solid and dot- 
ted curves correspond to solutions using the upper sign in 
Eq. (|37|l and dashed curves to solutions with the lower sign 
in Eq. (|37[) . Panel a) shows a high-frequency branch, the EM 
branch shifted ±/i'/4 compared to the plasma frequency at 
fco = 0, and a low-frequency branch. For H > 4/3\/3 ~ 0.77 
the pair branch merge with the upshifted EM branch, and 
there is an instability for waves with small wavenumbers; the 
growth rate indicated with the dotted curve for H = 1. Panels 
d)-f) shows the dispersion curves for finite amplitude (ao = 1) 
EM waves. 



exhibits an instability for small wavenumbers. We have 
plotted the grovsfth rate for the instability for if = 1 
in Fig. [TJ;. Direct numerical simulations of the Dirac- 
Maxvsrell system have confirmed this instability. It leads 
in the nonlinear stage to an interplay between the \E'i and 
^4 components of the spinor, and the excitation large 
amplitude (ao = 1~2) oscillatory EM fields. The disper- 
sion relation for a finite amplitude EM wave, shown in 
Fig. [TJi)-f), shows that the frequency is downshifted in 
the intense EM field, and that the quantum frequency 
shifts decrease compared to the linear cases. Equation 
P7p with the lower sign also yields a low-frequency wave, 
plotted in Figs.[T^)-{TJ:) for Wq ^ c^fcp ^ oj^g, from which 
we have the low-frequency dispersion relation 



Wo 






(38) 



which is a low-frequency spin-EM wave. We found that 
the low-frequency branch exists as a propagating wave 
only in the weakly relativistic regime, and disappears 
completely as a propagating wave for ao ^ 0.2 AH. We 
mention that ion dynamics can become important in 
the low-frequency range. We mention that ion dynam- 
ics can become important in the low-frequency range. 
For cold fiuid ions, Eq. (P^ is replaced by Wq 



^2l2 



fCn 



-{Xe+Xi)^Q where the ion susceptibility is Xi 



-w. 



pt 



H 



and ujpi is the ion plasma frequency, and w^j is added 
to the right-hand side of Eq. 



pt 



we retain ([38] 

r,2 



il.2 



for a;2^, 

« ^p^ « 



Wn < C- 



'fcn < 



For this case, 
ojpg, while for 



UJq < c 

frequency ion mode loq 
mass. 



we instead have the low- 
fik'^/2mi where rrii is the ion 



pe' 



2. The dipole field 





a„=0.01 


5 




1 


4.5 


■ 


4 


■ 


3.5 


■ 


3 


■ \ 


|„2.5 


■ \ - 


2 


\ - 


1.5 


'J 


1 


^-~^^^_: 


0.5 





5r 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 
0.5 




5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 
0.5 





FIG. 2: The cutoff frequency tJo at if = fco = as a func- 
tion of H for different values of oo = eAo/meC. The cutoff 
frequency is upshifted for the spin 'up' state (solid lines) and 
downshifted for spin 'down' (dashed lines). In the classical 
limit ii — )• 0, we have cjo — >■ <^pe/y7o for both spin states, 
where 70 — a/I + oJ. 



It is also possible to find simple expressions for the 
plasma susceptibility for the dipole case k^ ^ K = 
with arbitrary amplitude ^Oj which yields the nonlinear 
cutoff frequency of the EM wave. Here, inserting ^1*2 — 
^3=0 into (IT4|) - pT)) yields nontrivial solutions for fi = 

Equations p^ or 



(l/n)v/(TOeC2 - htOo/2Y + c2e2A2. 

(fT7| then yield the relation between ^1 and ^4, which are 
normalized such that '^\ + '^1 = uq. On the other hand, 
inserting ^^ = 1^4 = into (IT1])-((T7I) yields nontriv- 
ial solutions for fl = {l/h)^J{mec'^ + ?iwo/2)2 -|- c^e^A^. 
Equations ([T5|) or ([T6)) then yield the relation between ^2 
and ^3, which are normalized such that Vfl -I- ^3 = n^. 
The resulting susceptibility is Xe — Xe±, where 



U}„ 



Xe± 



'^oTi' 



and we have denoted 



7. - J(l T P^Y + ^'^° 



2mpC^ 



mr.c^ 



(39) 



(40) 



which, inserted into (|24l) with fco = 0, yields 



IV. EXTENSIONS TO MIXED STATES AND 
KINETIC MODELS 



Wo 



UJ„ 



7± 



(41) 



Here ujq can be seen as the effective plasma frequency in 
the presence of quantum spin effects and the EM field. 
We note that the spin effect contributes to a relative fre- 
quency shift of the order ±?iWpe/4"ZeC^ compared to the 
classical plasma frequency, while a large amplitude ra- 
diation field Aq leads to a frequency downshift, which 
resembles the effect of the relativistic electron mass in- 
crease in the classical plasma. 

The dependence of luq on H has been plotted in Fig. [5] 
for different values of ag = eAo/meC. We see that the 
frequency shifts increase linearly with H for H <^ 1, while 
the upshifted branch experiences a sharp rise at H fa 
0.77, which is the critical value of H where the upper 
branch looses its stability according to Fig. [T] We note 
that the relative quantum shift disappears both in the 
classical limit ?i — > and in the non-relativistic limit 
c — >■ oo, hence the relative shift TiLUpel^nieC^ is a combined 
quantum and relativistic collective effect. 



3. Vacuum case 

It is also interesting to consider the vacuum case 
Wo = cfco, originally considered by Volkov [36], where 
we have either D+ = or £)_ = in Eq. (IT51) . 
We consider here A as an external field, not influ- 
enced by the plasma, and calculate the plasma re- 
sponse. For £)_!_ = 0, the quasi-neutrality and current- 
neutrality conditions (^(1)) and ([^ give K — —ko/2 — 
—ui^jlc and 17 = — a;o/2 + ^^711^0^ /fi, and the solu- 
tions *i = [(l-|-7o)/27o]J'o, *2 = -(eAo/TOec) $0/270, 
^3 = -(e2A2/m2c2)^o/27o(l + 7o), and ^4 = 
(eAo/mec)^o/27o. If we instead use Z?_ = 0, then 
we obtain K — fco/2 = wo/2c, and il — ujq/2 + 
-fonieC^/h, and the solutions $1 = (eAo/mec)$o/27o, 
^2 = [(l + 7o)/27o]*o, *3 = (eAo/mec)$o/27o, and 

The resulting susceptibility for both cases is 



^^070' 



(42) 



which is identical to the case of the classical plasma where 
the relativistic gamma factor gives rise to nonlinear ef- 
fects, such as the self-induced transparency of the EM 
waves [30[. The above result was also obtained in a sim- 
plified model [4l| by using the Klein-Gordon equation for 
spinless particles. Hence, for the vacuum case, there is 
no difference in plasma response between the two spin 
states. 



In the above investigation, we have considered the ide- 
alized case where all electrons have a well-defined spin 
state. Therefore, these results can be seen as limiting 
cases of more complicated cases with an admixture of 
electrons with different spin. The simplest mixed state 
solution could be achieved by assuming that we have an 
admixture of the electrons with spin-up and spin-down 
states. If the electrons are distributed equally among the 
two spin states (but not among negative energy states), 
we would instead of ([M|) have the dispersion relation 



2 _ 2,2 _ 



{Xe++Xe-) 



^0, 



(43) 



where Xe+ and Xe- are the electric susceptibilities ob- 
tained by solving the two separate Dirac equations for 
the two spin states, each normalized according to ((20)) 
and fulfilling ([22]). For example, for the linear case we 
would use the susceptibilities Xe± in ([35]) in (|43|) to con- 
struct the dispersion relation 



(44) 
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which can be rewritten as 
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(45) 



For ujQ w ujpe and cfco ^ Wpe, there is a relative quan- 
tum upshift of the frequency of the order H'^ /&, which is 
extremely small for normal laboratory conditions where 
H '^ 1. On the other hand, for Wg <C c^k^ ■ 
have a low- frequency branch loq — hc^kQ/2uipemeC 
ing a similar procedure for the nonlinear cutoff frequency, 
the dispersion relation (|4T|) would then be replaced by 



';e> we 
■2 Us- 



^pe f I 1 



^ — + — =0, 

2 V7+ 7-/ 



(46) 



where 7± are given by (|40p . More complex, kinetic mod- 
els can be constructed by extending the dispersion rela- 
tion (|43l) to include not only the sum over the two spin 
states, but also sums (or integrals) over the excited states 
having different wavenumbers K. We like to mention 
that a kinetic model recently has been derived [43| to de- 
scribe wave propagation in a relativistic quantum plasma 
based on the Klein-Gordon equation and using a Wigner 
transform technique. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the nonlinear propaga- 
tion of light in dense matter by using a collective Dirac 
model coupled with the Maxwell equations, which in- 
cludes the nonlinear effects of the finite amplitude EM 



waves and the electron spin-1/2 effects. As an example, 
we have considered the nonlinear propagation of circu- 
larly polarized EM waves in a quantum plasma, and have 
studied the effects of different spin polarization, which in- 
troduces an up- or down shift of the EM waves, depend- 
ing on whether the plasma electron is in a spin 'up' or 
spin 'down' state. This relative frequency shift is of the 
order 10~^-10~^ for typical solid density or compressed 
density plasmas in the laboratory, but could be much 
larger in astrophysical settings (e.g. in the core of white 
dwarf stars), where the plasmonic energy density htdp,, is 
comparable to the electron rest energy. The spin related 
frequency shift could potentially be observed experimen- 
tally if a high-density plasma slab [40| with a definite 
spin state is irradiated with a laser beam with different 



polarization. In such a plasma it is also expected that 
linearly polarized laser light would perform Faraday ro- 
tation due to the different dispersive properties of the 
right-hand and left-hand polarized wave. For laser waves 
with frequencies between the cutoff frequencies for the 
right-hand and left-hand polarized wave, the plasma slab 
would work as a filter and only allow one of the polar- 
izations to propagate through the slab. Above a critical 
plasma number density, there is a density driven insta- 
bility, in which the pure spin up or spin down state is 
unstable under the generation of circularly polarized EM 
waves. Instabilities of this type could be important in 
the core of white dwarf stars, if the plasma has been spin 
polarized by a strong magnetic field. 



J 



Appendix A: The Dirac equation for circularly polarized EM waves 

For the case of circularly polarized EM waves of the form A = y4o(xcos6' — ysin6') where 9 — k^z — wot, we seek 
solutions of the Dirac equation ([T]) of the form ip — '^{d)exp{iKz — iftt), where we introduced the wavenumber K 
and frequency fl that are related to the momentum and energy of the electrons. This yields the relations 



-^ = - ( tJo-^ + «f^* ) exp(i/s:z - iVLt), 



(Al) 



and 



Vi/" = z ( fco h iK^ \ expiiKz - iflt). 

' d9 ' 



(A2) 



The Dirac equation ([T} takes the form 



ihiujQ h iil'^) — COL ■ 

do 
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m^c P^ = 0, 



(A3) 



which is a coupled system of 4 ordinary differential equations for ^1-^4. To put it in an explicit form, we evaluate 
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(A6) 



The Dirac equation thus takes the form 



[ - iH^o-7^ + i^) - mec^] *i + ich{ko— + iK)-93 - ceAoe'^*4 = 
[ - ih{uJo— + ifl) - nieC^] *2 - ich,{ko— + iK)'^4: - ce^oe"^^*3 = 

Lilly \A/\j 

[ - M^a-j^ + i^) + mec^] *3 + ich{ko— + iA')*i - ceAoe'^*2 = 

I" - ihiwo—- + iVl) + meC^l *4 - ich{ko—- + iK)^2 ~ ceAoe^'^^i = 
^ dO ■' dO 



(A7) 

(A8) 

(A9) 

(AlO) 



To ehmmate the e*^ and e~*^ phase factors, we assume *i = ^le'''/^, vj/j = ^2e~''^^'^, *3 = *3e*^/^, and *4 = 
^46"*^/^, where ^1-^4 are constants, to obtain 



[n{n + ^) - mec^j^i - c;i(X + ■^)*3 - ce^o*4 = 



2 



fco. 



[n(r! - y) - m,c^]^2 + ch{K - ^)*4 - ceAo*3 = 



Wo- 



ko. 



[h{n + y) + mec2] *3 - ch{K + ^)*i - ceylo*2 = 



W0^ 



^0. 



[h{n - y) + mec2]*4 + c;i(X - ^)*2 - ce^o^i = 0, 



(All) 
(A12) 
(A13) 
(A14) 



where J7 takes the role of an eigenvalue. The coefficient matrix is real and symmetric, hence fl is real and ^j can also 
be assumed real. 



Appendix B: Derivation of the electron current 



The X, y and 2 components of the current density 

je — —ecip^aip = —ec{ip''ax'4'^ + '>P cty'>Py + "0 "zV'^)j 
are obtained with the help of the expressions (See Appendix |^ for the relation between ipj and \I'j) 
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= ($1*4 + *3*2)(e*'' + e-''^) = 2(^1*4 + *3$2) cos(6'). 
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= ^(*l*4 + *3*2)(e''' - e-*^) = -2(*i*4 + ^3*2) sin(6l). 
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respectively, giving 



je = -2ec[(v&i*4 + vI/3V]/2)(xcos6'- ysin^) + (*i*3 - *2*4)z] 
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